Visualizing the Central Nervous System: Imaging Tools for Multiple Sclerosis and Neuromyelitis Optica Spectrum Disorders

Multiple sclerosis (MS) and neuromyelitis optica spectrum disorders (NMOSD) are autoimmune central nervous system conditions with increasing incidence and prevalence. While MS is the most frequent inflammatory CNS disorder in young adults, NMOSD is a rare disease, that is pathogenetically distinct from MS, and accounts for approximately 1% of demyelinating disorders, with the relative proportion within the demyelinating CNS diseases varying widely among different races and regions. Most immunomodulatory drugs used in MS are inefficacious or even harmful in NMOSD, emphasizing the need for a timely and accurate diagnosis and distinction from MS. Despite distinct immunopathology and differences in disease course and severity there might be considerable overlap in clinical and imaging findings, posing a diagnostic challenge for managing neurologists. Differential diagnosis is facilitated by positive serology for AQP4-antibodies (AQP4-ab) in NMOSD, but might be difficult in seronegative cases. Imaging of the brain, optic nerve, retina and spinal cord is of paramount importance when managing patients with autoimmune CNS conditions. Once a diagnosis has been established, imaging techniques are often deployed at regular intervals over the disease course as surrogate measures for disease activity and progression and to surveil treatment effects. While the application of some imaging modalities for monitoring of disease course was established decades ago in MS, the situation is unclear in NMOSD where work on longitudinal imaging findings and their association with clinical disability is scant. Moreover, as long-term disability is mostly attack-related in NMOSD and does not stem from insidious progression as in MS, regular follow-up imaging might not be useful in the absence of clinical events. However, with accumulating evidence for covert tissue alteration in NMOSD and with the advent of approved immunotherapies the role of imaging in the management of NMOSD may be reconsidered. By contrast, MS management still faces the challenge of implementing imaging techniques that are capable of monitoring progressive tissue loss in clinical trials and cohort studies into treatment algorithms for individual patients. This article reviews the current status of imaging research in MS and NMOSD with an emphasis on emerging modalities that have the potential to be implemented in clinical practice

Imaging markers of disability in aquaporin-4 immunoglobulin G seropositive neuromyelitis optica: a graph theory study

Neuromyelitis optica spectrum disorders lack imaging biomarkers associated with disease course and supporting prognosis. This complex and heterogeneous set of disorders affects many regions of the central nervous system, including the spinal cord and visual pathway. Here, we use graph theory-based multimodal network analysis to investigate hypothesis-free mixed networks and associations between clinical disease with neuroimaging markers in 40 aquaporin-4-immunoglobulin G antibody seropositive patients (age = 48.16 ± 14.3 years, female:male = 36:4) and 31 healthy controls (age = 45.92 ± 13.3 years, female:male = 24:7). Magnetic resonance imaging measures included total brain and deep grey matter volumes, cortical thickness and spinal cord atrophy. Optical coherence tomography measures of the retina and clinical measures comprised of clinical attack types and expanded disability status scale were also utilized. For multimodal network analysis, all measures were introduced as nodes and tested for directed connectivity from clinical attack types and disease duration to systematic imaging and clinical disability measures. Analysis of variance, with group interactions, gave weights and significance for each nodal association (hyperedges). Connectivity matrices from 80% and 95% F-distribution networks were analyzed and revealed the number of combined attack types and disease duration as the most connected nodes, directly affecting changes in several regions of the central nervous system. Subsequent multivariable regression models, including interaction effects with clinical parameters, identified associations between decreased nucleus accumbens (β = −0.85, P = 0.021) and caudate nucleus (β = −0.61, P = 0.011) volumes with higher combined attack type count and longer disease duration, respectively. We also confirmed previously reported associations between spinal cord atrophy with increased number of clinical myelitis attacks. Age was the most important factor associated with normalized brain volume, pallidum volume, cortical thickness and the expanded disability status scale score. The identified imaging biomarker candidates warrant further investigation in larger-scale studies. Graph theory-based multimodal networks allow for connectivity and interaction analysis, where this method may be applied in other complex heterogeneous disease investigations with different outcome measures

Uncovering convolutional neural network decisions for diagnosing multiple sclerosis on conventional MRI using layer-wise relevance propagation

Machine learning-based imaging diagnostics has recently reached or even superseded the level of clinical experts in several clinical domains. However, classification decisions of a trained machine learning system are typically non-transparent, a major hindrance for clinical integration, error tracking or knowledge discovery. In this study, we present a transparent deep learning framework relying on convolutional neural networks (CNNs) and layer-wise relevance propagation (LRP) for diagnosing multiple sclerosis (MS). MS is commonly diagnosed utilizing a combination of clinical presentation and conventional magnetic resonance imaging (MRI), specifically the occurrence and presentation of white matter lesions in T2-weighted images. We hypothesized that using LRP in a naive predictive model would enable us to uncover relevant image features that a trained CNN uses for decision-making. Since imaging markers in MS are well-established this would enable us to validate the respective CNN model. First, we pre-trained a CNN on MRI data from the Alzheimer's Disease Neuroimaging Initiative (n = 921), afterwards specializing the CNN to discriminate between MS patients and healthy controls (n = 147). Using LRP, we then produced a heatmap for each subject in the holdout set depicting the voxel-wise relevance for a particular classification decision. The resulting CNN model resulted in a balanced accuracy of 87.04% and an area under the curve of 96.08% in a receiver operating characteristic curve. The subsequent LRP visualization revealed that the CNN model focuses indeed on individual lesions, but also incorporates additional information such as lesion location, non-lesional white matter or gray matter areas such as the thalamus, which are established conventional and advanced MRI markers in MS. We conclude that LRP and the proposed framework have the capability to make diagnostic decisions of..

Comparison of probabilistic tractography and tract-based spatial statistics for assessing optic radiation damage in patients with autoimmune inflammatory disorders of the central nervous system

Background: Diffusion Tensor Imaging (DTI) can evaluate microstructural tissue damage in the optic radiation (OR) of patients with clinically isolated syndrome (CIS), early relapsing-remitting multiple sclerosis and neuromyelitis optica spectrum disorders (NMOSD). Different post-processing techniques, e.g. tract-based spatial statistics (TBSS) and probabilistic tractography, exist to quantify this damage. Objective: To evaluate the capacity of TBSS-based atlas region-of-interest (ROI) combination with 1) posterior thalamic radiation ROIs from the Johns Hopkins University atlas (JHU-TBSS), 2) Juelich Probabilistic ROIs (JUEL-TBSS) and tractography methods using 3) ConTrack (CON-PROB) and 4) constrained spherical deconvolution tractography (CSD-PROB) to detect OR damage in patients with a) NMOSD with prior ON (NMOSD-ON), b) CIS and early RRMS patients with ON (CIS/RRMS-ON) and c) CIS and early RRMS patients without prior ON (CIS/RRMS-NON) against healthy controls (HCs). Methods: Twenty-three NMOSD-ON, 18 CIS/RRMS-ON, 21 CIS/RRMS-NON, and 26 HCs underwent 3 T MRI. DTI data analysis was carried out using JUEL-TBSS, JHU-TBSS, CON-PROB and CSD-PROB. Optical coherence tomography (OCT) and visual acuity testing was performed in the majority of patients and HCs. Results: Absolute OR fractional anisotropy (FA) values differed between all methods but showed good correlation and agreement in Bland-Altman analysis. OR FA values between NMOSD and HC differed throughout the methodologies (p-values ranging from p < 0.0001 to 0.0043). ROC-analysis and effect size estimation revealed higher AUCs and R squared for CSD-PROB (AUC=0.812; R squared=0.282) and JHU-TBSS (AUC=0.756; R squared=0.262), compared to CON-PROB (AUC=0.742; R squared=0.179) and JUEL-TBSS (AUC=0.719; R squared=0.161). Differences between CIS/RRMS-NON and HC were only observable in CSD-PROB (AUC=0.796; R squared=0.094). No significant differences between CIS/RRMS-ON and HC were detected by any of the methods. Conclusions: All DTI post-processing techniques facilitated the detection of OR damage in patient groups with severe microstructural OR degradation. The comparison of distinct disease groups by use of different methods may lead to different - either false-positive or false-negative - results. Since different DTI post-processing approaches seem to provide complementary information on OR damage, application of distinct methods may depend on the relevant research question

Differences in Advanced Magnetic Resonance Imaging in MOG-IgG and AQP4-IgG Seropositive Neuromyelitis Optica Spectrum Disorders: A Comparative Study

Aims: To explore differences in advanced brain magnetic resonance imaging (MRI) characteristics between myelin oligodendrocyte (MOG) immunoglobulin (IgG) and aquaporin-4 (AQP4) IgG seropositive (+) neuromyelitis optica spectrum disorders (NMOSD). Methods: 33 AQP4-IgG and 18 MOG-IgG seropositive NMOSD patients and 61 healthy control (HC) subjects were included. All 112 participants were scanned with the same standardized MRI-protocol on a 3-Tesla MRI-scanner. Brain volume and diffusion tensor imaging (DTI) parameters were assessed. Results: MOG-IgG+ patients showed reduced parallel diffusivity within white matter tracts compared to HC whereas AQP4-IgG+ showed no significant brain parenchymal damage in DTI analysis. AQP4-IgG+ patients showed reduced whole brain volumes and reduced volumes of several deep gray matter structures compared to HC whereas MOG-IgG+ patients did not show reduced brain or deep gray matter volumes compared to HC. Conclusions: Microstructural brain parenchymal damage in MOG-IgG+ patients was more pronounced than in AQP4-IgG+ patients, compared with HC, whereas normalized brain volume reduction was more severe in AQP4-IgG+ patients. Longitudinal imaging studies are warranted to further investigate this trend in NMOSD. Our results suggest that MOG-IgG+ and AQP4-IgG+ NMOSD patients differ in cerebral MRI characteristics. Advanced MRI analysis did not help to differentiate between MOG-IgG+ and AQP4-IgG+ patients in our study

MRI Markers and Functional Performance in Patients With CIS and MS: A Cross-Sectional Study

Introduction: Brain atrophy is a widely accepted marker of disease severity with association to clinical disability in multiple sclerosis (MS). It is unclear to which extent this association reflects common age effects on both atrophy and function. Objective: To explore how functional performance in gait, upper extremities and cognition is associated with brain atrophy in patients with Clinically Isolated Syndrome (CIS) and relapsing-remitting MS (RRMS), controlling for effects of age and sex. Methods: In 27 patients with CIS, 59 with RRMS (EDSS <= 3) and 63 healthy controls (HC), 3T MRI were analyzed for T2 lesion count (T2C), volume (T2V) and brain volumes [normalized brain volume (NBV), gray matter volume (NGMV), white matter volume (NWMV), thalamic volume (NThaIV)]. Functional performance was measured with short maximum walking speed (SMSW speed), 9-hole peg test (9HPT) and symbol digit modalities test (SDMT). Linear regression models were created for functional variables with stepwise inclusion of age, sex and MR imaging markers. Results: CIS differed from HC only in T2C and T2V. RRMS differed from HC in NBV, NGMV and NThaIV, T2C and T2V, but not in NWMV. A strong association with age was seen in HC, CIS and RRMS groups for NBV (r = -0.5 to -0.6) and NGMV (r = -0.6 to -0.8). Associations with age were seen in HC and RRMS but not CIS for NThaIV (r = -0.3; r = -0.5), T2C (r(s) = 0.3; r(s) = 0.2) and T2V (r(s) = 0.3; r(s) = 0.3). No effect of age was seen on NWMV. Correlations of functional performance with age in RRMS were seen for SMSW speed, 9HPTand SDMT (r = -0.27 to -0.46). Regression analyses yielded significant models only in the RRMS group for 9HPT, SMSW speed and EDSS. These included NBV, NGMV, NThaIV, NWMV, logT2V, age and sex as predictors. NThalV was the only MRI variable predicting a functional measure (9HPT(r)) with a higher standardized beta than age and sex (R2 = 0.36, p < 1e-04). Conclusion: Thalamic atrophy was a stronger predictor of hand function (9HPT) in RRMS, than age and sex. This underlines the clinical relevance of thalamic atrophy and the relevance of hand function as a clinical marker even in mildly disabled patients

Functionally Relevant Maculopathy and Optic Atrophy in Spinocerebellar Ataxia Type 1

Background: Spinocerebellar ataxia type 1 (SCA-ATXN1) is an inherited progressive ataxia disorder characterized by an adult-onset cerebellar syndrome combined with nonataxia signs. Retinal or optic nerve affection are not systematically described. Objectives: To describe a retinal phenotype and its functional relevance in SCA-ATXN1. Methods: We applied optical coherence tomography (OCT) in 20 index cases with SCA-ATXN1 and 22 healthy controls (HCs), investigating qualitative changes and quantifying the peripapillary retinal nerve fiber layer (pRNFL) thickness and combined ganglion cell and inner plexiform layer (GCIP) volume as markers of optic atrophy and outer retinal layers as markers of maculopathy. Visual function was assessed by high- (HC-VA) and low-contrast visual acuity (LC-VA) and the Hardy-Rand-Rittler pseudoisochromatic test for color vision. Results: Five patients (25%) showed distinct maculopathies in the ellipsoid zone (EZ). Furthermore, pRNFL (P < 0.001) and GCIP (P = 0.002) were reduced in patients (pRNFL, 80.86 ± 9.49 μm; GCIP, 1.84 ± 0.16 mm3) compared with HCs (pRNFL, 97.02 ± 8.34 μm; GCIP, 1.98 ± 0.12 mm3). Outer macular layers were similar between groups, but reduced in patients with maculopathies. HC-VA (P = 0.002) and LC-VA (P < 0.001) were reduced in patients (HC-VA [logMAR]: 0.01 ± 010; LC-VA [logMAR]: 0.44 ± 0.16) compared with HCs (HC-VA [logMAR]: –0.12 ± 0.08; LC-VA [logMAR]: 0.25 ± 0.05). Color vision was abnormal in 2 patients with maculopathies. Conclusions: A distinct maculopathy, termed EZ disruption, as well as optic atrophy add to the known nonataxia features in SCA-ATXN1. Whereas optic atrophy may be understood as part of a widespread neurodegeneration, EZ disruption may be explained by effects of ataxin-1 gene or protein on photoreceptors. Our findings extend the spectrum of nonataxia signs in SCA-ATXN1 with potential relevance for diagnosis and monitoring

7 Tesla MRI of Balo's concentric sclerosis versus multiple sclerosis lesions

Background: Baló’s concentric sclerosis (BCS) is a rare condition characterized by concentrically layered white matter lesions. While its pathogenesis is unknown, hypoxia-induced tissue injury and chemotactic stimuli have been proposed as potential causes of BCS lesion formation. BCS has been suggested to be a variant of multiple sclerosis (MS). Here, we aimed to elucidate similarities and differences between BCS and MS by describing lesion morphology and localization in high-resolution 7 Tesla (7 T) magnetic resonance imaging (MRI) scans. Methods: Ten patients with Baló-type lesions underwent 7 T MRI, and 10 relapsing remitting MS patients served as controls. The 7 T MR imaging protocol included 3D T1-weighted (T1w) magnetization-prepared rapid gradient echo, 2D high spatial resolution T2*-weighted (T2*w) fast low-angle shot and susceptibility-weighted imaging. Results: Intralesional veins were visible in the center of all but one Baló-type lesion. Four Baló-type lesions displayed inhomogeneous intralesional T2*w signal intensities, which are suggestive of microhemorrhages or small ectatic venules. Eight of 10 BCS patients presented with 97 additional lesions, 36 of which (37%) had a central vein. Lesions involving the cortical gray matter and the U-fibers were not detected in BCS patients. Conclusion: Our findings support the hypothesis that BCS and MS share common pathogenetic mechanisms but patients present with different lesion phenotypes

Optic chiasm measurements may be useful markers of anterior optic pathway degeneration in neuromyelitis optica spectrum disorders

Objectives: We aimed to evaluate optic chiasm (OC) measures as potential imaging marker for anterior optic pathway damage assessment in the context of neuromyelitis optica spectrum disorders (NMOSD). Materials and method: This cross-sectional study included 39 patients exclusively with aquaporin 4-IgG seropositive NMOSD of which 25 patients had a history of optic neuritis (NMOSD-ON) and 37 age- and sex-matched healthy controls (HC). OC heights, width, and area were measured using standard 3D T1-weighted MRI. Sensitivity of these measures to detect neurodegeneration in the anterior optic pathway was assessed in receiver operating characteristics analyses. Correlation coefficients were used to assess associations with structural measures of the anterior optic pathway (optic nerve dimensions, retinal ganglion cell loss) and clinical measures (visual function and disease duration). Results: OC heights and area were significantly smaller in NMOSD-ON compared to HC (NMOSD-ON vs. HC p < 0.0001). An OC area smaller than 22.5 mm2 yielded a sensitivity of 0.92 and a specificity of 0.92 in separating chiasms of NMOSD-ON from HC. OC area correlated well with structural and clinical measures in NMOSD-ON: optic nerve diameter (r = 0.4, p = 0.047), peripapillary retinal nerve fiber layer thickness (r = 0.59, p = 0.003), global visual acuity (r = - 0.57, p = 0.013), and diseases duration (r = - 0.5, p = 0.012). Conclusion: Our results suggest that OC measures are promising and easily accessible imaging markers for the assessment of anterior optic pathway damage

Foveal changes in aquaporin‐4 antibody seropositive neuromyelitis optica spectrum disorder are independent of optic neuritis and not overtly progressive

Background and purpose: Foveal changes were reported in aquaporin-4 antibody (AQP4-Ab) seropositive neuromyelitis optica spectrum disorder (NMOSD) patients; however, it is unclear whether they are independent of optic neuritis (ON), stem from subclinical ON or crossover from ON in fellow eyes. Fovea morphometry and a statistical classification approach were used to investigate if foveal changes in NMOSD are independent of ON and progressive. Methods: This was a retrospective longitudinal study of 27 AQP4-IgG + NMOSD patients (49 eyes; 15 ON eyes and 34 eyes without a history of ON [NON eyes]), follow-up median (first and third quartile) 2.32 (1.33-3.28), and 38 healthy controls (HCs) (76 eyes), follow-up median (first and third quartile) 1.95 (1.83-2.54). The peripapillary retinal nerve fibre layer thickness and the volume of combined ganglion cell and inner plexiform layer as measures of neuroaxonal damage from ON were determined by optical coherence tomography. Nineteen foveal morphometry parameters were extracted from macular optical coherence tomography volume scans. Data were analysed using orthogonal partial least squares discriminant analysis and linear mixed effects models. Results: At baseline, foveal shape was significantly altered in ON eyes and NON eyes compared to HCs. Discriminatory analysis showed 81% accuracy distinguishing ON vs. HCs and 68% accuracy in NON vs. HCs. NON eyes were distinguished from HCs by foveal shape parameters indicating widening. Orthogonal partial least squares discriminant analysis discriminated ON vs. NON with 76% accuracy. In a follow-up of 2.4 (20.85) years, no significant time-dependent foveal changes were found. Conclusion: The parafoveal area is altered in AQP4-Ab seropositive NMOSD patients suggesting independent neuroaxonal damage from subclinical ON. Longer follow-ups are needed to confirm the stability of the parafoveal structure over time